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摘 要 : 在 宽带 


法 。 该 算法 首先 通 


情况 下 ， 


关键 词 : 


波束 形成 中 ， 若 干扰 信号 
化 。 为 了 解决 波束 形成 中 的 这 
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些 问题 ， 
过 EMP 算法 求 取 阻 塞 和 矩阵 ， 
空间 (CSM) 方法 、 协 方差 给 阵 重 构 求 取 合理 的 协 方差 矩阵 ; 最 后 进 
该 算法 能 
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研究 了 一 种 基于 协 方差 矩阵 重 构 和 特征 投影 预 处理 (EMP) 的 宽带 波 
村 干扰 相 消 预 处 理 阻 塞 掉 主办 干扰; 然后 通过 相 


对 接收 信号 进 和 


波束 形成 的 波 


从 主 鸭 方向 进入 ， 会 引起 主办 畸变 、 穷 狂 电 平 抬 高 ， 从 而 使 波束 性 能 严重 恶 


束 形 成 算 
干 信号 子 


提 行 波束 形成 。 在 主办 干扰 、 旁 办 干扰 同时 存在 的 
够 自 适 应 地 阻塞 主办 和 干扰、 抑制 劣 六 和 干扰， 解决 了 存在 主办 干扰 情况 下 宽带 
题 。 计 算 机 仿真 验证 了 该 算法 有 效 性 

宽带 波束 形成 ; 主因 干扰 ; 特征 投影 预 处 理 ; 协 方差 矩阵 重 构 


形 畸 变 问 


Broadband beamforming algorithm based on covariance matrix reconstruction and EMP 
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Abstract: Mainlobe interference can cause the distortion of the mainlobe and the increase of the sidelobe level in broadband 


beamforming. This paper developed a broadband beamforming algorithm based on EMP (eigen-projection processing) algorithm 


and covariance matrix reconstruction. This algorithm constructed a blocking matrix, preprocessed the echo signal. , then obtained 


an appropriate covariance matrix, finally formed beam. This algorithm could block the mainlobe interference and suppressed 


the sidelobe interference and solved the problem of waveform distortion in the presence of mainlobe interference and sidelobe 


interference. Simulation verify the effectiveness of the algorithm. 
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0 引言 
域 上 9。 近 些 年 来 许多 学 


E 带 波束 形成 算法 广泛 应 用 于 雷达 、 
者 针对 汰 兴 干扰 提出 了 很 多 抑制 劳 办 


通信 、 电 子 对 抗 等 领 


干扰 的 波束 形成 算法 。 文 献 [和 讨论 了 稀 疏 阵列 的 抑制 旁 久 干扰 


算法 , 文献 [5] 研 究 了 基于 干扰 噪声 


Capon 算法 提高 了 输出 信 干 品 比 且 算 法 具有 稳健 性 。 
成 中 若干 扰 信 号 从 主 洲 
EE 问 题 ;: a) 蔗 办 日 


钊 


具有 更 高 的 输 


波束 
现 两 个 玫 
-对 上 述 


阵 预 处 型 
的 波束 形成 。 


和 矩阵 Capon 波束 形成 使 阵列 
信和 干 品 比 ， 文 献 [6,7] 分 别 应 用 压缩 感知 和 稳健 


问题 提 


法 : a) 基 于 阻塞 矩 


方向 进入 ， 波 束 形成 就 会 出 
EE 平 抬 高 ; b) 主 波束 严重 变形 。 文献 [8] 
出 了 两 种 阻塞 主 欠 干扰 的 方 


(BMP ) 的 波束 
其 中 BMP 算法 存在 一 些 缺 陷 : 


成 ; b) 基 于 特 和 


投影 预 处 理 (EMP) 


旁 汶 较 高 » 并 且 对 


主办 干扰 的 估计 角度 要 求 严格 ;而 EMP 算法 也 存在 波峰 偏 移 
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男 ， 四 川 源 人 ， 


的 缺点 。 
进 ,解决 了 波峰 1 


文献 [9] 就 文献 [8]JEMP 算法 的 一 些 问题 作出 了 一 些 改 


局 移 的 问题 .对 于 主办 干扰 下 宽带 


波束 的 研究 


文献 [10] 给 出 了 一 种 处 理 方法 , 但 是 该 方法 旁 浴 电 平 很 高 , 这 将 
导致 虚 警 概率 的 急剧 上 升 。 
本 文 针 对 文献 [8] 中 EMP 算法 存在 的 波峰 


了 协 方差 矩阵 重 构 ， 并 且 将 


EMP 算法 应 用 到 了 宽 


1 ” 宽 囊 信和 号 模型 


假设 有 N 个 均匀 分 布 的 阵 元 。 信 号 


局 本 当 人 沼 
源 总 个 数 


有 相同 的 带 
宽 。 阵 列 输 


过 快速 离 


砚 忆 ， 噪声 为 高 共 
出 数据 疝 量 x(t) 如 下 : 


i 移 的 问题 进行 
带 信 号 上 。 


为 Dp 个 ， 它 们 具 


X(t) = [xi (2), ,xn (OT 


散 侍 里 叶 变 换 (DFT) 分 


i 


折 白 噪声 且 与 源 信号 具有 相同 的 带 


(1) 


} 解 为 / 子 频带 。 在 中 心 


河南 灵 宝 人 ， 


副教授 ， 博 士 ， 主 要 研究 方向 为 电子 对 抗 、 高 速 实时 信号 处 


录用 稿 


X(fi)=AF): SF) tNF) j=12%, (2) 

其 中 : 为 第 /个 子 频带 信号 的 中 心 频率 ; X(f) 为 N x 1 子 带 接收 
号 向 量 ; 8( 廊 ) 为 信号 在 第 /个 子 频带 的 谱 分 量 ，N(f) 为 Nx 1 
为 噪声 在 户 处 的 谱 向 量 ; 4()= 
a(gn, 亡 )] 为 第 万 个 子 频带 内 信和 号 的 方向 


e-j2r(d/DGw-Dsinei]T， 


,D，D 为 信号 源 


向 量 ， 
[a(01,fi), a(02,f),, 
矩阵 ，a(6i fi) = [1,e-j27(4/1) sin Bi... 
9; 为 第 i 个 信号 相对 于 阵 元 的 入 射 角 ， 


2 ” 子 频 带 干 扰 相 消 预 处 理 


假设 空间 中 有 P 个 旁 辩 方向 干扰 信号 ， 一 个 主流 方向 干扰 
信号 ， 则 第 1/ 个子 频带 只 包含 干扰 和 噪声 的 协 方差 矩阵 为 Rx(j)， 
特征 值 分 解 如 下 : 
Rx( 万 ) = Ed =UAUS + UmnAmUR (3) 
其 中 : XGi = 12,...,N) 是 协 方差 矩阵 Rx (fi) 的 特征 值 ， 且 A4 > 
?和 2 二 …>Mp>Mpri>…hvi 全 为 与 特征 值 丰 对 应 的 特征 向 量 ， 
…,Up+1 张 成 干扰 子 空间 ，Up42， Up+3,…,UN 张 成 噪声 
子 空 间 ， 而 Us = [wuz Up4i), A, = 
diag (M1,h2,°%* ,Ap41), Un = [Upi42, Up43, UN An = 
diag(4p+2,4p+3,…,AN)。 在 实际 中 ，Rx( 记 ) 可 以 用 实际 接 
收 信 号 的 协 方差 矩阵 访 x( 太 ) 代 替 ， 


让 2 


U1, U2,* 


Rx(f) = EK X(f) XS) (4) 


kK 表示 训练 快 拍 数 。 
则 第 /个 信道 的 阻塞 矩阵 BUDD 可 以 用 式 〈3) 表示 为 
B(fi)=IT— un unm) lum® (5) 
其 中 : wu 表 示 主 瓣 干扰 所 对 应 的 特征 向 量 , 并 且 I 是 N x N 的 单位 
阵 。wum 由 式 《6) 确定 : 


Jp (wm a(05,f))| = maxlo (wal0,f))| 《6) 


其 中 : p(qa1,q2) 表 示 问 量 aj， az 的 相关 系数 被 定义 为 


0 0 (7) 


llaillllazll 

其 中 ，a(8。 亡 ) 表 示 期 望 信号 在 方向 9, 频率 万 上 的 导向 向 量 。 

对 中 心 频率 为 / 子 频带 接收 信号 干扰 相 消 预 处 理 后 信号 为 
Yb(f;) = B(f)X(f) (8) 

: X(f) 可 以 看 成 多 个 空间 信号 合成 ， 即 
J = a(0s,fi)s(0s,f;) + Pe a(0i,f;)s(0i, fi) + N(f) 
(9) 

其 中 :s(6s, 万 ) 表 示 期 望 信号 在 /个 信道 的 复 包 络 ,s(gp+l 万 ) 

表示 主 瓣 干扰 在 /个 信道 的 复 包 络 ，s(6u 户 )(0 = 1,2,…,P) 表 示 

旁 状 干扰 的 复 包 络 。 式 (10) 为 经 过 干扰 相 消 预 处 理 后 的 信号 

即 


p(a1, Q2z) 全 


Y5(f) = 》[a(eb 方 ) -umGnanon-Tamnaa(gb 方 sb 有 


i=1 


+[a(05, fi) — um (um um) um a(0s, f;)|s(0s, f;) 


毛 寄 达 ， 等 ,一 种 基于 协 方差 矩阵 重 构 和 EME 的 宽带 人 


+[a(Op+1,f;) — Umum um) um a(Op+1f;)]s(Op+1,f;) 
+N(f) (10) 
其 中 : WP) = B(f;)N() 为 预 处 理 后 的 噪声 。 由 于 ww 为 主 拆 干 
扰 对 应 的 特征 向 量 可 知 
Um!a(0,fi)~ 0 i= 1,2,.…P (11) 
Um(Um! um) lumia(Op+1,fi) ~ a(0p+1f/) (12) 
可 知 预 处 理 后 子 频 带 的 信号 向 量 为 


Ys(f) = a(0f))s(05,f) + > a(0uf))s(01,f) + NA) 
i=1 


(13) 

可 以 看 出 经 过 子 频 带 干扰 相 消 预 处 理 后 , 第 j 个 子 频带 阵列 

接收 数据 中 的 主办 干扰 成 分 被 阻塞 , 劳 浴 干扰 成 分 被 保留 下 来 ， 
方便 下 一 部 分 波束 形成 。 


3 ”聚焦 处 理 


>s 间 ) 算法 求解 聚焦 
VvV( 有 为 
9s)A"(fi,0s) 的 


本 算法 采用 文献 [11] 中 的 RSS (旋转 子 空 
矩阵 。 在 RSS 算法 中 T(j) =V( 有 040f)， 其 
A(fo,0s)47(f,0s) 的 左 奇异 向 量 ; U( 记 ) 为 4(ho， 
右 奇 异 向 量 , ,为 期 望 信号 的 方向 。 其 中 : 


H 
A(f,0:) = er ee] (14) 


A 0) = [eer (sme, .oI2rc-D(E) smo (15) 
聚焦 后 第 /个 信道 最 终 的 输出 为 
Y( 让 =T(P)(B(CP)X(CP)) (16) 
此 时 的 信号 协 方差 矩阵 为 
Ry(f) =Y(J)Y2 (万 ) (17) 


4 ， 协 方差 矩阵 重 构 


首先 如 知 采 用 文献 [12] 中 经 典 的 SMI《〈 和 采样 矩阵 求 逆 ) 算 
法 求 取 权 向 量 ， 第 /个 子 频 带 自 适 应 权 回 量 为 


OD a 
1 a(Osf) HRYy(F) a(0sf)) 


设 4 = 1/a(0s, 有 j)”Ry(fj) a(0s, 有 j)， 经 过 对 协 方差 和 
行 特征 值 分 解 ， 式 〈18) 可 以 写 为 


N 
和 
wi = 4) ,Tua(0sf) 
i=1 


(18) 


| 
种 


车 进 


过 -|o(e。 万 ) 一 2 uu a(g,, 万]| (19) 


并 求 其 自 适 应 方向 图 得 到 
Fa(9,f) 加 wj"a(0,f;) 


a {Fa(0,f) — Der ie [al0s, f) | Fi(0,#))} 


(20) 
Fa 人 6, 万) 为 第 /个 子 频带 中 的 自 适 应 方向 图 :Fa(6, 万) 为 静 
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.hiNaXn 站] 


作 期 
录用 稿 i 达 ， | 皮 束 形 la 


态 方向 图 ; Fi(9, 及) 为 特征 波束 , 且 Fo(9, fi) = a(0s, fj)?a(0, ff)， 


Fi(0,f;) = ui*a(0, fi), Fi(0,fi)( = 12,…,P+1) 为 P +1 个 干 
扰 对 应 的 特征 波束 ，Fi(9, 万) 在 第 i 个 干扰 信号 方向 上 的 增益 最 
大 , 在 其 他 干扰 方向 增益 为 零 ; 4min 为 协 方差 矩阵 Rr (万 ) 的 最 小 
特征 值 。 

当天 为 干扰 信号 (= 1,2,… 
大 于 Amin， 因 此 可 以 得 到 


Ni 
1 


min > 1 


,P 二 1) 对 应 的 特征 值 时 ，Xi 要 


i = 1,., (21) 


可 以 看 出 特征 向 量 % 在 它 对 应 干扰 方向 上 形成 了 零 陷 。 
若 将 N 设 为 噪声 对 应 的 特征 值 时 人 me 
近 ， 因 此 


=P+2,P+3,. 


Ni 


一 (22) 


所 以 假如 把 主办 干扰 特征 值 改 为 噪声 特征 值 ， 这 样 求 出 的 
权 疝 量 受 到 主流 干扰 的 影响 ， 从 而 避免 波峰 偏 移 的 问题 。 在 本 
文 算法 中 ， 如 式 (23) 所 示 ， 主 办 干扰 对 应 的 特征 值 ji 设 为 噪 


声 特征 值 的 均值 各， 从 而 将 主办 干扰 特征 值 降 为 背景 噪声 的 水 
平 ， 主 办 干扰 特征 值 对 权 向 量 影响 就 不 存在 了 。 
An = 02 一 Tt (23) 
a 羽 为 噪声 对 应 的 特征 值 会 产 
生 扰 动 ， 这 会 使 天 线 方向 图 旁 办 抬 高 ， 影 响 波 束 形 成 性 能 ， 所 
以 为 了 消除 噪声 抖动 ， 将 协 方差 矩阵 噪声 特征 值 修改 为 
Tn a (24) 
可 以 得 出 第 /个 子 频 带 协 方差 矩阵 可 以 表示 为 
Ry(f) = U402 (25) 
其 中 : 4 为 重 构 后 的 特征 值 矩阵 ， 且 
A= diag(41,%, dm Ap+1, G2 7T2) (26) 
本 算法 中 特征 向 量 并 没有 发 生 改 变 ， 经 ; 过 襄 焦 处 理 后 协 方 
差 矩阵 为 
Ry =37/17(f)Ry TS) (27) 
根据 MVDR (最 小 方差 无 畸变 响应 ) 算法 ， 求 得 宽带 波束 
成 的 权 向 量 关 
W(fo) = LRy a0s, fo) (28) 
其 中 :f= 1/a(9,, 用)?UA-1UHa(9,, 万 ) 为 一 标量 。 可 以 看 出 在 


求 取 权 癌 量 时 ， 主 办 干扰 特征 值 的 影响 已 经 消除 ， 因 此 权 癌 量 
不 会 使 天 线 方向 图 产生 波峰 偏 移 。 A 为 对 角 阵 , 所 以 对 其 求 逆 仅 
仅 只 是 对 角 元 素 的 倒数 。 因 此 不 用 进行 复杂 的 矩阵 求 逆 运 算 ， 
运算 量 比 矩阵 求 逆 降 低 了 许多 。 阵 列 最 终 频 域 输出 为 Z(f) = 
[WO) YR), ,WY ), .WY(f)] 时 域 为 z(t) = 
IDFT[Z(f)]。 

本 文 算法 信号 的 处 理 流程 如 图 1 所 示 。 


5 仿真 实验 及 


本 次 实 


号 ， 本 次 仿真 信 


两 个 为 旁 瓣 干扰 。 


旁 锥 干扰 来 波 方向 分 别 为 40`、 
通过 DFT 分 为 / = 128 个 子 频 带 。 图 2~4 均 为 三 
维 图 在 最 高 频率 fy = 2.1GHz、 


阵列 接收 的 信号 


1.9GHz 处 的 截面 图 。 
图 2 为 未 经 过 


阵 重 构 求 
取 权 值 


图 1 


本 文 算法 信号 处 理 流程 


结论 


验 采 用 一 个 阵 元 数 为 16, 阵 元 间距 为 工作 频带 最 高 
频率 对 应 的 半 波 长 。 假 设 期 望 信号 


、 和 干扰 信号 均 为 线性 调频 信 


号 中 心 频率 有 = 2GHz， 带 宽 B = 200MHz。 了 噪声 为 
高 斯 白 噪 声 , 信 品 比 为 0dB。 有 三 个 干扰 信和 号， 


个 为 主办 干扰 、 
主 泥 干 扰 来 波 方 向 为 1, 干 噪 比 为 40dB。 两 个 
一 30", 干 品 比 分 别 为 20dB、20dB。 


! 心 频率 f= 2GHz、 最 低频 率 fi = 


过 干扰 相 消 预 处 理 直 接 采用 MVDR 的 仿真 图 。 


可 以 看 出 天 线 方向 图 在 主 辩 干扰 方向 2 附近 出 现 了 很 深 的 零 陷 ， 
这 严重 影响 了 波峰 的 位 置 而 且 旁 辩 很 高 。 


图 


图 3 所 示 为 经 过 


02 
角度 0 ( 度 ) 


2 MVDR 算法 形成 的 天 线 方向 图 


EMP 干扰 相 消 预 处 理 ,但 未 进行 协 方差 矩 


阵 重 构 ， 最 后 采 月 


日 MVDR 的 仿真 图 。 可 以 看 出 EMP 算法 构建 


的 阻塞 矩阵 阻塞 掉 了 主 准 干扰， 天 线 方向 图 并 未 在 主办 干扰 方 
向 形成 零 陷 ， 但 是 波峰 位 置 发 生 了 2 的 偏 移 。 
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1 增益 (dB) 


-50 02 50 
角度 0 ( 度 ) 


EMP 算法 


图 3 采用 成 的 天 线 方向 图 


图 4 为 基于 协 方差 矩阵 重 构 和 EMP 的 宽带 波束 形成 仿真 
图 .可 以 看 出 本 文 算法 成 功 解决 了 EMP 算法 的 波峰 偏 移 问题 。 
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图 4 本 文 算 法 的 天 线 方向 


表 1 给 出 了 三 种 算法 波束 指向 、 主 办 宽度 、 旁 闪电 平 的 具 
体 数 值 对 比 。 


表 1 各 种 算法 性 能 的 对 比 


MVDR EMP 本文 算 法 
波束 指 据 —3.7° 2 0.1° 
fo 处 主办 宽度 (-3dB) 5.7° 9° 9.2° 
fi 处 主因 宽度 (-3dB) 6.0° 9.6" 9.7° 
全 处 主办 宽度 (-3dB) 5.5° 8.6" 8.8- 
和 劳 激 电 乎 值 —3.13dB —12.9dB  —13.9dB 


实验 证 明 , 对 于 经 典 MVDR 算法 并 不 能 阻塞 主办 干扰 ， 所 
以 从 仿真 图 中 可 以 看 出 主办 干扰 零 陷 严 重 影响 期 望 信号 的 接收 
且 旁 闪电 平 很 高 波形 严重 恶化 。EMP 算法 可 以 阻塞 主办 和 干扰， 
但 是 可 以 看 出 波峰 位 置 发 生 了 偏 移 。 本 文 算法 通过 协 方差 矩阵 
重 构 解决 了 波峰 偏 移 的 问题 ， 而 且 只 需要 特征 值 的 分 解 就 能 求 
出 子 频 带 阻 塞 和 矩阵 B0) 和 权 值 wh)， 并 且 BG)、W0W) 可 以 并 行 
运算 ， 提 高 算法 的 执行 效率 。 经 典 MVDR 算法 复杂 度 为 0(] * 
K3),， 本 文 算法 和 EMP 算法 复杂 度 为 o(2 *J* K3)。 本 文 算法 解 
决 了 经 典 MVDR 算法 的 波形 畸变 问题 ， 相 比 EMP 算法 在 不 改 
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变 算 法 复杂 度 的 情况 下 解决 了 波峰 偏 移 问题 。 


6 ”结束 语 


本 文采 用 基于 协 方差 矩阵 重 构 和 EMP 的 算法 ， 对 存在 主 
办 干扰 信号 情况 下 的 宽带 波束 形成 算法 进行 了 研究 。 本 文 算法 
分 为 三 个 部 分 处 理 ，a) 各 个 子 频带 求 取 阻塞 矩阵 阻 寒 掉 主办 干 
扰 ; b) 采 取经 典 RSS 算法 对 接收 信号 进行 聚焦 处 理 ，c) 对 聚焦 
后 的 协 方差 矩阵 进行 重 构 ， 和 矩阵 重 构 解决 了 波峰 偶 移 问题 。 经 
过 计算 机 仿真 可 以 看 出 ， 本 文 算法 可 以 处 理 在 主办 干扰 情况 下 
的 宽 这 波束 形成 ,尤其 是 当主 流 干 扰 严 重 影 响 天 线 方向 图 时 候 ， 
可 以 做 到 天 线 方向 图 的 保 形 。 
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